Jet modification in 200 AGeV Au-Au collisions 
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The computation of hard processes in hadronic collisions is a major success of perturbative Quan- 
tum Chromodynamics (pQCD). The environment of heavy-ion collisions offers the opportunity to 
embed such hard processes into a soft medium which is created simultaneously and study the 
medium-induced modifications. On the level of single high transverse momentum (Pt) hadrons, a 
reduction in yield, the so-called quenching is observed. However, on the level of jets, the energy- 
momentum flux carried by hadrons is conserved, i.e. the effect of the medium is a redistribution of 
energy and momentum, and statements about quenching of jets can only be made for specific cuts 
used to identify the jet. In this paper, we present a baseline computation for the expected modifica- 
tion of the longitudinal momentum distribution in jets by the medium created in 200 AGeV Au-Au 
collisions given a realistic set of experimental cuts used to identify jets in a heavy- iron environment. 
For this purpose, we use a Monte-Carlo code developed to simulate in-medium shower evolution 
following a hard process and a 3-d hydrodynamical evolution to simulate the soft medium. 



PACS numbers: 25.75.-q,25.75.Gz 



I. INTRODUCTION 



Jet quenching, i.e. the energy loss of hard partons cre- 
ated in the first moments of a heavy ion collision due to 
interactions with the surrounding soft medium has long 
been regarded a promising tool to study properties of 
the soft medium U 0, [3, |j, 0, @1 • The basic idea is to 
study the changes induced by the medium to a hard pro- 
cess which is well-known from p-p collisions. A number 
of observables is available for this purpose, among them 
suppression in single inclusive hard hadron spectra Raa 
, the suppression of back- to-back correlations @, or 
single hadron suppression as a function of the emission 
angle with the reaction plane • 

Single hadron observables and back-to-back correlations 
are well described in detailed model calculations using 
the concept of energy loss [lj]) E3, > i-e- under the as- 
sumption that the process can be described by medium- 
induced radiation which shifts the leading parton energy 
by an energy AE, followed by a fragmentation process us- 
ing vacuum fragmentation with the shifted energy. How- 
ever, there are also calculations for these observables in 
which the evolution of the in-medium parton distribu- 
tion is followed in an analytic way [IH, [l5|, [l6j . Recently, 
also Monte Carlo (MC) codes for in-medium shower evo- 
lution have become available [lT], EH- In such compu- 
tations, energy is not simply lost but redistributed in a 
characteristic way. 

Jet observables (as opposed to single hadron measure- 
ments) allow in principle to probe this redistribution of 
energy through the interaction with the medium in detail. 
In this paper, we aim to make a benchmark calculation 
for the modification of the longitudinal momentum dis- 
tribution of jets in 200 AGeV central Au-Au collisions 



at RHIC, based on t he j et finding strategy used by the 
STAR collaboration fl9j |. This observable in particular 
is expected to be sensitive to energy redistribution to 
increased hadron production at low Pt due to medium- 
enhanced branchings [20| . 



II. MODELLING JETS IN MEDIUM 

We describe jet production in several steps. First we 
compute the production of hard partons in leading order 
(LO) pQCD. These primary partons subsequently evolve 
as a parton shower and we assume that it is this partonic 
evolution which is modified by the medium. Eventually 
the parton shower hadronizes, and we assume that this 
process takes place sufficiently far outside the medium to 
be treated as in vacuum. The hadronized parton shower 
constitutes a jet in our model and is analyzed using the 
experimental jet finding strategy. A detailed description 
of the model is found in [l8j . here we just summarize the 
main steps. 

The production of two hard back to back partons k, I with 
momentum px in a p-p or A-A collision in LO pQCD is 
given by 



da 



AB^kl+X 



dp T dy 1 dy 2 



= 22xifi/A(xi,Q 2 )x2f ] /B{x2,Q 2 ) 



da* 



df 
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(1) 

where A and B stand for the colliding objects (protons 
or nuclei) and 2/1(2) is the rapidity of parton k(l). The 
distribution function of a parton type i in A at a mo- 
mentum fraction X\ and a factorization scale Q ~ px is 
fi/A(xi,Q 2 ). The distribution functions are different for 
the free protons [U, [22j and protons in nuclei [H, [24[ • 
The fractional momenta of the colliding partons i, j are 
given by x 1<2 = ^ (exp[±j/i] + exp[±y 2 ]). 
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Expressions for the pQCD subprocesses — (s, i, u) as 
a function of the parton Mandelstam variables s, t and u 
can be found e.g. in [25J. To account for various effects, 
including higher order pQCD radiation, transverse mo- 
tion of partons in the nucleon (nuclear) wave function 
and effectively also the fact that hadronization is not a 
collinear process, we fold into the distribution an intrinsic 
transverse momentum fcy with a Gaussian distribution of 
width 1.6 GeV. This momentum vector points into a ran- 
dom direction in the transverse plane and modifies the 
transverse momenta px t , Pt 2 °f the outgoing partons in 
creating a momentum imbalance between them according 
to p Tl +Pt 2 = k T . 

The probability density P(xo,yo) for finding a hard ver- 
tex in an A-A collision at the transverse position ro = 
(xo,Uo) and impact parameter b is given by the product 
of the nuclear profile functions as 

P / \ ?Mro + b/2)?Mro - b/2) 

where the thickness function is given in terms of 
Woods-Saxon the nuclear density p A (r,z) as T^(r) = 
/ dzpA^.z). For the present study, we evaluate Eq. @ 
at b = corresponding to central collisions. 
Inclusive production of a parton flavour / at rapidity yj 
and momentum pt is found from Eq. {T]) by integrating 
over either y\ or ?/2 and summing over appropriate com- 
binations of partons. We use this distribution probed at 
midrapidity as an input for the next step, the evolution 
of a shower. 

In vacuum, we simulate using the PYSHOW algorithm 
[H which is part of PYTHIA [27j]. To take into account 
medium effects, we assume a medium which does not ab- 
sorb momentum by recoil of its constituents, but rather 
increases the virtuality of partons propagating through it 
and thus inducing additional radiation. Such a medium 
can be characterized by a transport coefficient q which 
represents the increase in virtuality AQ 2 per unit path- 
length of a parton traversing the medium. 
Since the pQCD shower evolution takes place in momen- 
tum space, in order to incorporate this type of medium 
modification, we need to make a link to the spacetime 
evolution of the shower. We assume that the formation 
time of a shower parton with virtuality Q is developed 
on the timescale l/Q, i.e. the lifetime of a virtual par- 
ton with virtuality Qb coming from a parent parton with 
virtuality Q a is in the rest frame of the original hard col- 
lision (the rest frame of the medium may be different by 
a flow boost as the medium may not be static) given by 

t, — E b _ E b 

~ Qb or 

The time r° at which a parton a is produced in a branch- 
ing can be determine by summing the lifetimes of all an- 
cestors. Thus, during its lifetime, the parton virtuality 
is increased by the amount 

&Ql = / d(q(C) (3) 



where C is integrated along the spacetime path of the 
parton through the medium and q(Q = q(T,r,(f>,r) s ) de- 
scribes the spacetime variation of the transport coeffi- 
cient where q is specified at each set of coordinates proper 
time t, radius r, angle 4> and spacetime rapidity rj s . In 
the following, we will use a hydrodynamical evolution 
model of the medium [28| for this dependence assuming 
that q scales as 

= J K'-2-e 3/4 (C)(coshp-sinlipcosa) (4) 

with the energy density e and the local flow rapidity p 
with angle a between flow and parton trajectory. If the 
parton is a gluon, the virtuality transfer from the medium 
is increased by the ratio of gluon to quark Casimir color 
factors, 3/| = 2.25. 

The additional virtuality transfer to the parton in Eq. ([3]) 
leads to an increased branching rate and thus to both 
to an increase of the angular spread of the shower with 
respect to the original hard parton trajectory and a re- 
duction in hard parton production due to increased pro- 
duction of soft partons [la] . Note that Eq. (J3]) needs 
to evaluated for every parton path through the medium 
from the hard vertex, Eq. |2]). However, in 18] it was 
found that to good approximation it is sufficient to sort 
paths according to AQ 2 ot , the integrated virtuality along 
the eikonal path of the original hard parton through the 
medium. This simplifies the computation considerably. 
In a last step, the Lund string fragmentation scheme [291 ] 
is used to hadronize the shower. Note that, as discussed 
in [la ], while the assumption that hadronization can be 
assumed to be unmodified by the medium is expected 
to hold for high Pt light hadrons, it may break down 
for heavy hadrons produced at low momentum, thus this 
end of the distribution is computed less reliably. 

III. JET IDENTIFICATION AND 
MEDIUM-INDUCED MODIFICATION 

Experimentally, jets are identified using calorimetric 
measurements in addition to a series of cuts to elimi- 
nate background contributions. Thus, for a given parton 
energy E p in the simulation, a detector does not always 
find a jet with the same energy. This introduces a bias on 
any properties of identified jets as compared to all jets. In 
the following, we study this problem simulating as closely 
as possible the cuts used by the STAR collaboration to 
find jets [l!|, although we stress that the formalism is 
completely general and could be used for any other set 
of cuts. In this section, we focus on dN/dEE , the distri- 
bution of the detected jet energy E, given an underlying 
parton energy E p . 

For an ideal detector, one would expect this distribution 
to be 6(E — E p ) in the absence of intrinsic kr and a Gaus- 
sian centered around E p when taking into account inrin- 
sic kx- However, the STAR calorimetric measurement is 
only sensitive to charged hadrons, ttq and 7, i.e. when- 
ever the hadronization yields a different neutral hadron, 
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FIG. 1: (Color online) Distribution of the energy with which a quark jet of 30 GeV energy is detected dN/dE, given different 
cuts in hadron momenta or cone width. Left panel: without medium modification. Right panel: medium modified with 
AQtot = 15 GeV 2 (see text), roughly corresponding to a path through the whole medium. 



its energy does not contribute to the reconstructed jet 
energy. This induces a long tail towards low energies. 
Furthermore, in order to separate jets from the substan- 
tial background in a heavy-ion environment, the following 
additional requirements are imposed: (1) a seed hadron 
of 7.5 GeV or more is required to be in the shower 
as starting point for the jet finding algorithm (2) only 
hadrons with Pt > 2 GeV are taken into account and (3) 
an angular cut with angle with respect to the jet axis of 
R < 0.4 is imposed. The effect of these cuts on dN/dEE p 
is shown in Fig. [1] for a quark jet in vacuum (left panel) 
and for a jet with AQ 2 ot — 15 GeV 2 (right panel, cor- 
responding to a path traversing the whole medium), all 
for E p — 30 GeV. ft is clearly seen how the medium acts 
to reduce the energy of the jet seen by the detector by 
reducing the probability that a seed hadron is found, by 
transporting energy out of the cone region R < 0.4 and 
by increasing hadron production below the Pt cut. 
If we wish to study the medium suppression of jets ob- 
served in e.g. 25-30 GeV or 30-35 GeV momentum bins, 
we have to average the results of dN/dEE p over different 
E p and parton type with the weight factors determined 
by Eq. |T]) and in addition average over all different tra- 
jectories from the vertex of origin Eq. |2|) through the 
medium and then study which fraction of jets is found 
in a momentum bin, given the jet finding cuts, as com- 
pared to the situation in vacuum. In this way, a jet sup- 
pression factor R?a\ can be defined. For this particu- 
lar model of energy redistribution and the imposed cuts, 
R j f A is 0.43 ± 0.04 for the 25-30 GeV momentum bin and 
0.46±0.05 for the 30-35 GeV bin (statistical errors only). 
In order to study the change in the momentum distri- 
bution along the jet axis, we focus on dN/dt; where 
£ = \uE/Pt is the logarithm of the jet energy divided by 
the energy of a hadron in the jet. When comparing this 
quantity for a single parton with known energy with and 
without the medium, a depletion at low £ corresponding 
to the reduction of high Pt hadrons and an enhancement 
at high £ corresponding to low Pt hadron production is 



expected [H H HH . 

However, the fact that a jet needs to pass the cuts to 
be identified experimentally introduces a strong bias on 
dN/d£. Events with multiple hadron production be- 
low the Pt cut are strongly suppressed as compared to 
events with multiple hadron production just above the 
cut. Likewise, copious low Pt hadron production reduces 
the chance to find a hard seed hadron. 



— full calculation 

— no cuts 




FIG. 2: (Color online)The ratio of dN/d£ in central 200 AGeV 
Au-Au collisions over the vacuum value calculated for two 
different situations: Black only for jets which are found with 
an energy of 25-30 GeV within a cone of R < 0.4 and for 
hadrons with Pt > 2 GeV, but where the distribution was 
extracted without cut in Pt and for a wider cone R < 0.7. 
Red for fixed parton energy of 30 GeV assuming that all jets 
are detected without cuts. 

In Fig. [21 we present two different scenarios to illustrate 
this effect. First, we make a calculation close to the ex- 
perimental situation in which the jet needs to pass the 
cuts discussed in the previous section and in which then, 
again in analogy with the experimental procedure, dN / <i£ 
is extracted for these jets for all Pt in an increased cone 
with R < 0.7. In the second scenario, we assume an ideal 
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detector with no cuts imposed which allows to identify 
each produced jet regardless of its modification. 
In both scenarios, we average over all possible paths 
through the medium and compute the ratio of the dN / dt; 
in medium and dN/dt; in vacuum. It is clearly seen how 
the trigger bias changes the expectations for larger £. Ef- 
fectively, it prevents the observation of any enhancement 
in the region £ > 2 where the unbiased calculation shows 
a strong effect. 

On the other hand, the two scenarios are fairly similar 
towards £ = 0. In this region, the jet is effectively a 
single hadron carrying all the parton momentum. Such 
a jet is not influenced by any of the cuts and hence not 
biased. The fact that the value reached in the ratio is 
close to Raa for single hadron spectra is not an accident, 
as in this limit the suppression of single hard hadrons is 
probed. 

IV. DISCUSSION 

We have presented a benchmark calculation for the mod- 
ification of jet properties, here the longitudinal momen- 
tum distribution of hadrons inside the jet, due to a 
medium as created in heavy-ion collisions. We find in 
principle sizeable modifications, characteristic for the 



pattern of energy redistribution in the medium, i.e. a 
strong depletion of high Pt hadrons in the jet and an 
enhancement of low Pt hadron production. However, in 
a more realistic calculation close to the experimental jet 
finding strategy, a strong trigger bias prevents one from 
observing the enhancement. The depletion of high Pt 
hadrons in the jet however should still be visible (and 
agree with the single high Pt hadron suppression). 
While details of the calculation depend on the model used 
to simulate the medium effect, it should be stressed that 
the analysis framework is rather general and could be 
applied to any other model of in-medium jet modifica- 
tion (such as e.g. JEWEL [l7j]). Likewise, the change 
of event structure induced by the trigger bias appears to 
be generic rather than a peculiarity of the model used 
here. Finally, a comparison with data would clearly be 
valuable to study in more detail at which point model 
assumptions like hadronization outside the medium even 
for low Pt hadrons break down. 
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